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ABSTRACT 


The  unsteady  response  of  two  oceanic  density  fronts  to 
local  atmospheric  forcing,  using  combinations  of  wind  stress 
and  surface  heat  flux,  is  investigated  with  an  embedded 
mixed  layer-general  circulation  model.  The  adjustment  of 
the  frontal  structure  is  dependent  upon  the  wind  stress 
direction  and  whether  there  is  surface  heating  or  cooling. 
In  cases  of  an  applied  wind  stress  alone  where  denser  water 
is  transported  toward  less  dense  water,  the  frontal  struc¬ 
ture  diffuses,  the  mixed  layer  depth  deepens,  and  cross- 
frontal  mixing  occurs.  In  cases  where  less  dense  water  is 
transported  toward  denser  water,  the  frontal  structure  is 
preserved,  mixed  layer  depth  is  preserved  and  cross-frontal 
mixing  is  minimized.  The  addition  of  surface  heating  shal¬ 
lows  the  mixed  layer  and  inhibits  vertical  mixing.  Inertial 
oscillations  are  observed  ir.  the  across-front  velocity 
field. 
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INTROD  OCTIOH 


A.  EACKGBCDSC  AND  PURPOSE 

Curing  recent  years, there  have  been  a  number  of  scien¬ 
tific  studies  that  describe  the  nature  and  existence  of 
oceanic  fronts.  Cromwell  and  Reid  (1956)  identified  a  frcnt 
by  an  abrupt  horizontal  density  change  at  the  sea  surface. 
Roden  (1976)  srated  that  a  front  should  be  considered  to 
exist  wherever  any  oceanic  state  variable  (density,  tempera¬ 
ture  cr  salinity)  reaches  a  relative  maximum  value  in  its 
horizontal  gradient. 

The  motivation  for  studying  fronts  and  ocean  prediction 
has  teen  stressed  by  Roden  (1976)  ,  Elsberry  and  Garwood 
(1979),  and  Niiler  (1982).  Elsberry  and  Garwood  (1979)  spe¬ 
cifically  point  out  antisubmarine  warfare  applications, 
fisheries  management  and  a  resource  for  climate  research 
data,  as  viable  and  pertinent  reasons  for  further  study. 
Sany  important  physical  processes  which  are  accompanied  by 
attendant  air-sea  interactions  occur  in  frontal  regions. 
Latent  and  sensible  heat  transfer,  evaporative  and  precipi- 
tative  processes,  barcclinic  currents,  salt  fluxes,  mass  and 
momentum  transport,  inertial-internal  wave  activity  and  mix¬ 
ing  all  can  occur.  Oceanic  frontal  areas  are  areas  of  sound 
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speed  changes  and  often  are  areas  of  increased  biological 
activity. 

This  thesis  examines  the  t ime-depandeat  response  of  a 
numerically-  simulated  oceanic  front  to  local  atmospheric 
forcing.  The  model  used  in  this  work  incorporates  the  forc¬ 
ings  of  wind  stress  and  surface  hearing  with  a  20-level 
primitive  equation  ocean  circulation  model  with  embedded 
mixed  layer  dynamics. 

B.  CBSERV ATICNAL  STUDIES 

Cromwell  and  Eeid  (1956)  described  a  simple  case  of  an 
upper  layer  of  warm,  less  dense  water  overlying  a  layer  cf 
cool,  denser  water  in  the  ccean.  The  front  in  this  case  is 
a  density  front  formed  by  the  horizontal  temperature  gradi¬ 
ent.  Observations  shewed  surface  convergence  along  the 
front,  sinking  motion  in  the  immediate  vicinity  of  the  frent 
and  divergence  at  lower  depths.  They  observed  that  the  wind 
affected  the  front  by  mixing  the  upper  layer,  reducing  the 
horizontal  surface  temperature  gradient,  and  enhancing  the 
temperature  jump  at  the  base  of  the  mixed  layer. 

Voorhis  and  Hersey  (1964)  observed  a  wintertime  thermal 
front  in  the  Sargasso  Sea.  They  found  an  average  tempera¬ 
ture  gradient  of  ■  deg  C/10  km,  and  a  maximum  in  one 
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instance  cf  1.5  deg  C/5  km.  This  was  accompanied  by  a  60 
cm/s  jet  flowing  to  the  east  along  the  edge  of  the  front. 
The  depth  of  the  front  corresponded  to  the  mixed  layer 
depth.  The  water  on  either  side  of  the  front  was  found  to 
be  well-mixed. 

Katz  (1069)  continued  the  work  of  Voorhis  and  Hersey  and 
derived  temperature-salinity  correlations  for  a  front 
observed  in  the  Sargasso  Sea  in  May.  He  found  the  front  to 
be  a  separaticn  cf  two  distinct  water  masses,  with  a  tran¬ 
sition  zcne  cf  only  a  few  meters  in  the  vertical.  He 
deduced  that  water  transport  along  the  front  maintained  the 
sharp  transition.  Although  not  observed,  Katz  presumed  mix¬ 
ing  must  occur  at  the  interface,  but  the  along-front  current 
prevented  any  fcuild-up  of  well-mixed  waters  in  that  area. 
He  also  hypothesized  that  a  vertical  component  to  the  cur¬ 
rent  must  exist  at  the  interface.  Thus,  the  interface  would 
be  preserved  by  both  horizontal  and  vertical  motions,  and 
the  vertical  velocity  would  account  for  the  observed  accumu¬ 
lation  of  well-mixed  waters  below  the  front. 

Eang  (1973)  studied  a  front  at  the  southern  end  of  the 
Benguela  Current.  This  was  an  intense  front  with  a  horizon¬ 
tal  temperature  gradient  cf  1.6  deg  C/km  and  a  temperature 
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change  of  8-10  dag  C  across  it.  There  was  also  local 
upwelling  occurring  due  tc  an  offshore  transport.  As  might 
be  expected  from  Ekman  theory,  Bang  found  the  front  tc 
intensify  under  southerly  winds  (Ekmaa  transport  being  to 
the  left  of  the  wind  in  the  Southern  Hemisphere).  Under  a 
northwest  cale,  the  front  rapidly  weakened.  It  reintensi¬ 
fied  as  the  wind  steadied  from  the  south. 

Rcden  has  made  significant  and  extensive  contributicns 
to  the  literature  on  oceanic  fronts,  especially  on  the 
large-scale  fronts  in  the  Pacific  Ocsan.  In  a  series  of 
papers,  he  has  observed  and  examined  the  Pacific  subarctic 
front  (Roden  1975,  1977),  the  subtropical  front  (Roden  1974, 
1975,  1980),  the  doldrum  front  (Roden  1974,  1975),  and  the 
subarctic-subtropical  transition  zone  (Roden  1971,  1972). 
In  the  1971  paper,  Rcden  concentrated  on  the  subarctic-sub- 
tropical  transition  zone  and  found  that  the  wind  stress 
played  an  important  rcle  in  determining  the  location  of  heat 
and  salt  flux  divergence  zones.  Rcden  (1972)  focused  on  the 
temperature  and  salinity  fronts  located  at  the  boundaries  of 
the  transition  zone;  the  subarctic  front  to  the  north  and 
the  subtropical  front  (which  is  also  a  density  front)  to  the 
south.  The  origins  cf  bcth  fronts  were  found  to  be  related 
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to  the  wind  stress,  geostrophic  flow  fields  and  heat  and 
salt  flux  divergence.  The  doldrum  front  was  included  in  a 
study  of  the  subtropical  front  (Roden  1974).  Further  work 
on  all  three  fronts  was  conducted  by  Roden  (1975)  with  ref¬ 
erence  to  the  wind  and  to  the  energy  flux  fields.  Re 
observed  than  the  depth  of  the  nixed  layer  changed  dramati¬ 
cally  across  the  subarctic  front  in  winter  and  spring;  from 
100  m  on  the  north  side  to  300  m  on  the  south  side.  The 
subtropical  front  is  highly  dependent  on  the  wind  and  energy 
flux  fields  at  the  sea  surface,  as  the  front  occurs  in  an 
area  of  Eknan  transport  convergence.  The  doldrum  front  is  a 
front  due  to  a  salinity  gradient  only,  with  a  baroclinic 
eastward- flowing  current  at  the  surface  and  a  faster  west- 
ward-flowing  undercurrent. 

For  the  subarctic  front,  Roden  (1977)  found  atmospheric 
forcing  (that  is,  wind  stress,  radiative  heat  flux,  precipi¬ 
tation  and  turbulent  energy  fluxes)  to  be  the  primary  cause 
of  frontal  movement,  intensification  and  decay.  Roden  has 
shown  repeatedly  that  the  roles  of  wind  stress,  Ekman  trans¬ 
port  and  heating  are  inextricable  from  the  dynamics  of  fron- 
togenesis.  Roden  (1980)  again  examined  the  subtropical 
frontal  zone  (within  which  are  located  a  number  of 
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individual  fronts)  and  showed  that  Ekman  transport  led  to  a 
concentration  of  the  temperature  and  salinity  gradients  in 
the  upper  ocean.  Thus,  frontogenesis  was  due  to  horizontal 
convergence  and  confluence  of  the  flow  field. 

While  Koden's  work  has  dealt  primarily  with  frontogene¬ 
sis  ard  in  explaining  the  dynamics  of  several  identified 
North  Pacific  fronts,  his  findings  are  pertinent  to  this 
thesis,  i.e.,  local  atmospheric  forcing  plays  a  major  role 
in  the  dynamics  of  oceanic  fronts. 

C.  MCDEL  STUDIES 

Garvine  has  concentrated  on  a  modelling  approach  to 
fronts.  In  cne  of  his  earlier  papers  (1974),  he  investi¬ 
gated  the  dynamics  of  small-scale  fronts,  though  without 
regard  to  atmospheric  forcing.  Such  fronts  propagate  intc 
the  ambient  water  via  the  unbalanced  horizontal  pressure 
gradient.  Because  of  minimal  cross-frontal  mixing,  he  found 
strong  convergence  and  sinking  at  the  front.  He  later 
expanded  this  model  to  include  the  effects  of  wind  stress 
(Garvine  1979a,  1979b)  .  However,  because  of  the  need  to 
maintain  steady  state  in  the  model,  he  prescribed  wind 
stress  to  be  of  secondary  importance  relative  to  the  hori¬ 
zontal  pressure  gradient.  Buoyancy  and  energy  budgets  were 
added  later  to  the  model  (Garvine  1980). 
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Kao  (1980)  investigated  the  Gulf  Stream  as  a  large-scale 
density  front  in  the  upper  ocean.  The  heart  of  his  work  was 
performing  a  scale  analysis,  which  revealed  three  length 
scales:  an  inertial  or  deformation  length  scale,  a  buoyancy 
length  scale,  and  a  diffusive  length  scale.  Kao's  conclu¬ 
sion  was  that  the  front  was  maintained  by  the  cross-Gulf 
Stream  circulation,  but  his  model  included  no  wind  stress  or 
heating. 

Cushman-Rcisin' s  (198  1)  model  was  based  on  Roden's  work 
on  large-scale  frontogenesis,  and  it  does  have  some  similar¬ 
ities  to  the  model  used  in  this  thesis  because  it  includes 
mixed  layer  dynamics.  An  important  difference  is  in  the 
scales.  Cushman-Roisin  investigates  large-scale  and  long¬ 
term  spin-up  of  a  front  in  response  to  a  wind  stress  curl. 
The  model  in  this  thesis  is  used  to  investigate  local 
response  of  a  pre-existing  front  to  local  winds  and  heating. 
Surface  wind  stress  is  a  primary  means  of  forcing  in  bcth 
models,  but  Cushman-Roisin  applies  a  strong  wind  stress  curl 
in  the  fern  cf  a  wind  field  that  changes  direction  at  the 
latitude  of  the  zonally-oriented  front.  Be,  on  the  other 
hand,  spin  up  a  horizontally  uniform  wind  stress  over  six 
hours,  after  which  it  maintains  a  constant  value.  The 
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results  obtained  by  Cushaan-Roisin  depend  auch  upon  the  curl 
of  the  wind  stress.  Cushaan-Roisin  finds  the  nixed  layer 
depth  is  at  a  aisimui  in  the  Biddle  of  the  front,  while 
Ekman  downwelling  is  at  a  laxiaua.  This  is  probably  due  to 
a  lack  of  any  wind  stress  at  the  front  itself,  although  the 
curl  is  large. 

Caaerlengc  (1982)  studied  the  large-scale  response  of 
upper  ocean  frcnts,  specifically  the  Pacific  subarctic 
front.  His  aodel  has  a  variable-spaced  horizontal  grid  with 
highest  resolution  in  the  frontal  area,  vertical  mixing  was 
not  included  in  the  model.  Assuming  a  pre-existing  front, 
Caaerlengc  applied  different  wind  stresses,  and  found  that 
except  in  the  case  of  the  passage  of  a  strong  cyclone,  the 
effects  of  wind  forcing  were  limited  to  the  upper  150  m. 
One  case  of  a  uniform  wind  stress  did  show  Ekmar.  transport 
and  a  horizontal  displacement  of  the  front. 
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II.  DE  SCBI F1I0N  2E  IBS  PROBLEM 


A.  GENERAL 

The  response  of  a  pre-existing  oceanic  front  to  local 
atnospheric  forcing  (wind  stress  and  surface  heating)  over  a 
tine  period  of  two  days  is  examined.  An  embedded  mixed  lay¬ 
er-ocean  circulation  model  is  used,  and  two  different  fron¬ 
tal  structures  are  treated.  The  first  numerically-inserted 
front,  hereafter  referred  to  as  Front  1,  is  somewhat  of  an 
artifact  in  that  the  isotherms  become  horizontal  at  a  rela¬ 
tively  shallow  depth  (see  Figs.  1,  2,  and  3  for  the  initial 
conditions) ,  so  that  Front  1  is  manifested  only  within  the 
upper  60  m.  The  temperature  field  is  vertically  stratified 
belcw  approximately  70  m.  The  front  is  about  10  km  wide. 
The  maximum  horizontal  temperature  gradient  at  the  surface 
is  1.5  dag  C/1.6  km.  The  mixed  layer  depth  is  a  numerical¬ 
ly-inserted  approximation  to  what  the  actual  mixed  layer 
depth  would  be.  A  geostrcphic  along-front  velocity  is  cal¬ 
culated  frcm  the  initial  temperature  field.  It  has  a  maxi¬ 
mum  speed  cf  32  cm/s  with  its  axis  located  at  the  surface  on 
the  warm  side  cf  the  front.  The  across-front  geostrophic 
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Front  1  is  a  viable  representation  of  a 


velocity  is  zero, 
smaller  mescscale  upper  layer  front,  and  the  model  results 
will  be  important  in  this  context.  The  second  type  of 
front,  hereafter  referred  to  as  Front  2,  is  a  simulation  of 
a  more  typical  front,  similar  to  those  reported  by  acden 
(1980).  The  horizontal  temperature  gradient  at  the  surface 
is  1.5  deg  C/3.2  km.  The  front  is  not  only  evident  at  the 
surface,  but  also  its  effect  on  the  temperature  field 
extends  throughout  the  thermocline  of  the  model  (see  Figs. 
4,  5,  and  6  for  the  initial  conditions) .  The  mixed  layer 

depth  again  is  a  numerically- ins erted  approximation  based  on 
the  shape  cf  the  temperature  field.  It  is  felt  that  adjust¬ 
ments  which  occur  shortly  after  the  start  of  integration  of 
the  model  will  compensate  for  any  errors  in  the  initial 
field.  A  gecstrophic  along-front  velocity  is  calculated 
from  the  initial  temperature  field.  Its  maximum  value  is 
100  cm/s  and  the  axis  lies  at  the  surface  on  the  warm  side 
of  the  front.  The  acrcss-front  gecstrophic  velocity  is 
zero.  With  both  fronts  oriented  in  a  right-hand  coordinate 
system,  z  increases  upwards,  y  increases  to  the  right,  and  x 
increases  out  cf  the  page. 
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B.  HYPOTHESES 

It  is  exacted  that  the  direction  of  the  wind  relative 
to  the  orientation  of  the  front  plays  an  inportant  role  in 
the  near-frcntal  circulation  and  density  structure.  Ekman 
transport  across  the  front  should  cause  upwelling  and  down- 
welling  features  near  the  surface  mani festation  of  the 
front.  Jchannessen,  et  al  (1977)  observed  that  the  thermo- 
cline  through  the  Maltese  front  shallowed  in  certain  areas, 
which  indicated  upwelling,  and  also  that  the  thermocline 
spatially  deepened  after  isotherms  surfaced.  In  Johannes- 
sen's  case,  this  was  attributed  to  the  presence  of  an  eddy. 
If  the  water  Basses  had  been  of  similar  types,  it  would  have 
been  indicative  cf  downwelling. 

With  regard  tc  one- dimensional  nixing,  Niiler's  (1975) 
mixed  layer  scdel  used  an  impulsive  wind  stress  and  no  heat¬ 
ing,  and  showed  that  initially  the  mixed  layer  deepened  rap¬ 
idly.  Strong  inertial  actions  can  be  generated  which  car 
lead  to  a  large  increase  in  available  turbulent  energy  and 
thus  to  a  rapid  (within  the  first  half  pendulum-day)  deepen¬ 
ing  cf  the  mixed  layer.  flellor  and  Durbin  (1975)  also  used 
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an  impulsive  wind  stress  of  2  c»  /s  (approximately  2 
dynes/cm^)  and  no  heating.  They  showed  that  the  mixed  layer 
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Pigure  1.  Front  1  Initial  T?npe  rature.  Only  *he  mddle 
section  of  the  grid  and  the  upper  175  n  are 
shewn.  Contour  interval  is  0.5  deg  C. 
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Figure  3.  Frcnt  1  Initial  0-Geo strophic  Velocity,  contour 
interval  is  8  ca/s. 
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Figure  6. 


Freni  2  Initial  0-Geo strophic  Velocity,  Contour 
interval  is  20  ca/s.  Dashed  lines  indicate 
negative  values. 


deepened  and  that  the  temperature  jump  also  increased 
slightly.  In  a  two-dimensional  model,  De  Szoeke  (1980) 
examined  the  effects  of  wind  stress  only  in  his  version  of 
the  Niiler  (1975)  model  and  found  that  in  the  case  of 
constant  wind  stress,  the  mixed  layer  deepened  rapidly  in 
the  first  half  pendulum-day  (t*Pi/f),  and  then  slowed.  An 
additional  effect  of  the  wind  stress  is  the  possible  augmen¬ 


tation  of  a  guasi-ge cstrcphic  frontal  shear  (in 
front  direction)  which  could  increase  vertical 
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A  surface  heating -coo ling  cycle  can  lead  to  stratifica¬ 
tion  under  light  to  moderate  winds  which  will  tend  to  modu¬ 
late  any  aixing  tendencies  brought  about  by  wind-generated 
turbulence.  Niiler  (1975)  showed  that  surface  heating  halts 
the  deepening  cf  the  mixed  layer.  Me  expect  the  same  result 
in  our  model.  Hellor  and  Durbin  (1975)  investigated  a  case 
of  a  combined  constant  wind  stress  (2  dynes/cm^)  and  a 
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constant  beating/cooling  value  (H  =  ±.01  cm-K/s) .  For  "ha 
heating  case  they  found  a  shallowing  and  leveling  off  of  the 
depth  of  the  mixed  layer  and  an  ever-increasing  temperature 
jump  value  at  the  base  of  the  nixed  layer.  The  same  heating 
condition  with  a  lower  wind  stress  (1  dyne/cm^)  resulted  in 
a  shallower  nixed  layer  depth  and  the  establishment  of  a 
double  thermocline,  a  remnant  of  the  initial  nixed  layer  and 
the  formation  of  a  new  shallower  layer  lua  to  the  influence 
of  heating. 

It  is  expected  that  inertial  currents  across  the  front 
nay  play  an  important  role  in  the  across-f rontal  transport 
of  heat,  mass,  and  momentum.  Phasing  of  inertial  wind- 
driven  currents  in  conjunction  with  vertical  mixing  may 
result  in  lateral  entrainment  or  "capture”  of  buoyant  (warn) 
water  by  the  adjacent  cooler  water  during  horizontal  oscil¬ 
lations  across  the  front  plane. 

C.  METHOD 

To  test  the  above  hypotheses  and  to  establish  the  basic 
dependence  of  frontal  response  to  local  atmospheric  forcing, 
four  combinations  of  surface  boundary  conditions  are  inves¬ 
tigated.  These  are  listed  in  Table  I. 
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TABLE  I 

Atmospheric  Forcing  used  in  Model  Buns 

lx  fdyr.es/ca  )  ^  y  (dynes/cm  )  Q(deq  C-cm/s) 

Case  I  +1.0  0  0 

Case  II  -1-0  0  0 

Case  III  +1.0  0  -.004 

Case  17  -1.0  0  -.004 

In  the  cases  above,  Tx  is  the  along-front  (x-direction)  com¬ 
ponent  of  wind  stress,  is  the  cross-front  (y-directicn) 

component  of  wind  stress,  and  heat  flux  <2  is  applied  to  the 
ocean  uniformly  at  the  surface  (heat  flux  is  positive 
upwards,  so  the  negative  value  indicates  warming  of  the 
ocean  from  the  air)  .  All  four  cases  are  examined  for  both 
Front  1  and  Front  2. 

Cases  I  and  II  examine  the  effects  of  wind  stress  paral¬ 
lel  to  the  front  only,  which  will  eventually'  generate  a 
steady  Ekman  transport  perpendicular  to  the  front.  The  wind 
stress  of  1  dyne/cm2  corresponds  to  a  wind  speed  of  7-8  m/s. 

In  all  cases  the  wind  stress  is  spun-up  from  0  at  t=0  to  its 

2 

maximum  absolute  value  (1.0  dynes/cm  )  by  hour  6.  From  hour 
6  to  hour  48  it  remains  constant.  Such  a  short  spin-up  time 
is  net  unreasonable,  as  winds  can  change  on  a  shorter  time 
scale,  for  example,  when  an  atmospheric  cold  front  passes. 
Cases  III  and  IV  impose  a  constant  surface  heating  on  each 
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of  the  twc  preceding  cases.  As  indicated  by  the  earlier 
one-dimensional  mixed  layer  models,  we  expect  that  the  beam¬ 
ing  value  will  net  be  overwhelmed  by  the  wind  stress,  and 
that  the  effects  of  both  dynamic  and  thermodynamic  forcing 
will  be  evident. 

The  initial  conditions  for  Front  1  are  shown  in  Figs.  1, 
2,  and  3.  Warmer,  less  dense  water  resides  on  the  left,  and 
cooler,  denser  water  resides  on  the  right.  The  initial  con¬ 
ditions  for  the  velocity  are  determined  by  the  associated 
initial  density  field,  which  is  balanced  by  a  geostrophic 
along-front  current  directed  out  of  the  figure  (u  positive) 
(see  Fig.  3)  .  The  mixed  layer  depth  varies  from  a  uniform 
depth  in  the  far  field  tc  a  shallowing  through  the  front. 
The  model  is  run  for  a  40  hour  time  period  for  each  of  the 
specified  atmospheric  forcings  shown  in  Table  I. 

In  Case  I  we  expect  tc  see  water  from  the  cold  (right) 
side  transported  to  the  warm  (left)  side,  upwelling  on  the 
cold  side  at  the  front  and  dcwnwelling  due  to  convective  ar.d 
turbulent  mixing  cn  the  warm  side.  The  mixed  layer  depth 
should  deepen,  and  the  front  should  diffuse.  In  Case  II  we 
expect  the  warm  water  to  be  transported  on  the  surface  to 
the  cold  side,  but  convective  mixing  to  be  absent  due  to  the 
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stable  condition  of  the  less  dense  water  overriding  the  den¬ 
ser  water.  Opveiling  should  occur  on  the  warm  side  at  the 
front,  and  only  minor  dovnvelling  should  occur  or.  the  cold 
side  due  to  the  effect  cf  wind  stirring  alone.  In  both 
Cases  III  and  IV  the  heating  should  have  a  shallowing  influ¬ 
ence  on  the  mixed  layer  depth  tending  to  offset  deepening 
and  ether  effects  caused  fcy  turbulent  mixing. 

The  initial  conditions  for  Front  2  are  shown  in  Figs.  4, 
5,  and  6.  In  this  configuration,  the  warmer,  lass  dense 
water  resides  cn  the  right  and  the  cooler,  denser  water 
resides  on  the  left.  The  accompanying  along-front  geos- 
trcphic  current  is  directed  into  the  page  (u  negative)  with 
its  maximum  value  at  the  surface.  The  mixed  layer  depth 
follows  the  contours  of  the  upper  chermociine,  shallowing 
across  the  front.  The  model  is  initiated  and  integrated  in 
time  as  it  was  for  Front  1. 

In  Case  I  (wind  stress  positive,  directed  out  of  the 
figure) ,  we  expect  the  Ekaan  transport  to  be  toward  the 
left.  The  results  for  this  case  should  be  analagous  to  the 
expected  results  for  Case  II  applied  to  Front  1  as  previ¬ 
ously  noted.  Case  II  should  be  compared  to  Case  I  for  Front 
1  and  Cases  III  and  IV  should  be  compared  to  Cases  IV  and 
III  respectively  cf  Front  1. 
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The  inertial  period  fcr  our  model  i3 


20  hours 


at  the  latitude  cf  36  deg  N.  Considering  that  in  this  acdel 
wind  stress  dees  not  attain  its  steady  maximum  value  until  6 
hours  have  elapsed,  we  should  find  a  maximum  in  inertial 
motion  somewhat  later  than  half  an  inertial  period  after 
time  zero,  at  perhaps  t=  1 5  hours  or  so. 

D.  DESCRIEHCN  CF  THE  HODEL 

Tha  primary  documentation  for  the  model  is  in  Adamec,  et 
al  (1981).  In  that  paper,  the  model  is  cast  in  radial  coor¬ 
dinates  to  investigate  the  oceanic  response  to  a  stationary 
ra dially-symmetric  hurricane.  In  our  case,  the  model  is  in 
right-hand  Cartesian  coordinates  ar.d  it  is  two-dimensional 
(y  and  z) .  It  is  believed  that  this  model  is  particularly 
well-suited  tc  the  study  cf  upper  ocean  fronts  because  it 
incorporates  the  Garwood  (1977)  bulk  mixed  layer  model  into 
the  Haney  (1980)  primitive  equation  model.  The  ocean  is 
assumed  hydrcstatic  and  incompressible,  where  density  is  a 
function  of  temperature  alone.  The  Coriolis  parameter  is 
constant  and  there  are  no  fluxes  of  mass,  momentum  or  heat 
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normal  to  any  boundary  except  at  the  surface,  where  atmos¬ 
pheric  forcing  is  applied.  The  horizontal  grid  spacing  is 
125  m  and  the  vertical  grid  spacing  is  variable  over  20  lay¬ 
ers  from  6  a  at  the  top  layer  to  64  m  at  the  bottom  layer. 
The  domain  size  in  y  and  z  is  10  0  1cm  by  500  m,  respectively. 
Hith  the  frcnt  pre-existing  as  an  initial  condition,  the 
initial  u-ccmponent  of  velocity  is  in  geostrophic  balance, 
and  the  initial  v-component  is  initially  set  to  zero. 

Drawing  upcn  Adamec,  et  al  (1981),  highlights  of  the 
model  are  outlined  below.  Table  II  lists  the  pertinent 
model  variables  and  constants  used  in  the  governing 
equations. 

The  governing  equations  in  their  two* dimensional  form  are  as 
follows: 


U  -  -ft  -ft  * 

tv  -  &(<!•«•) 

(D 

ft  -  -ft  -  -ft  - 

*»  -  *ft  * 

-x(v.w., 

(2) 

ft  -  -ft  -  -ft  * 

^0  -  *(-*■) 

(3) 

(4) 
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TABLE  11 

flcdel  Constants  and  Variables 


At 

am 

Kip 

s 

t 

g 

fo 

To 


coefficient  of  thermal  expansion 
hcrizcntal  eijdy  conductivity 
coefficient  for  heat 
horizontal  eddy  diffusion 

coefficient  for  momentum 
vertical  eddy  conductivity 
coefficient  for  heat 
vertical  eddy  diffusion 

coefficient  for  aoaentun 
latitude 

Ccriclis  parameter 


2.0E-5  /degC 

2.5E+5  cm^/s 

2.5E+5  cra^/s 

0.5  cm^/s 

0.5  ca^/s 
36  N 

8.  55E-5/S 


gravity 

reference  density  of  seawater 


981.0  cn/s2 
1.0276  g/cm^ 


reference  temperature 


5  degC 


At  acdel  timestep 

J  horizontal  grid  array  size 

K  vertical  grid  array  size 

z=-D  basin  depth 

y  basin  width 

u  x  (along-front)  component  of 

velocity 

v  y  (across-front)  component  of 

velocity 

T  temperature 

h, H  mixed  layer  depth 


30  s 
800 
20 

500  a 
100  km 


0  "  "  If  "  (5) 

e  -  7  -  V  }  (*) 


The  boundary  conditions  are: 


u  *  0 


h 

May  " 


0  «  y  »  0 


100  km 


(7) 


w  *  0  @  z  *  0  ,  -D 


(8) 
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The  vertical  turbulent  fluxes  are  parameterized  in  two 
different  Banners  depending  upon  where  in  the  water  column 
they  are  located.  Eelow  the  mixed  layer,  the  vertical 
fluxes  are  parameterized  by  eddy  viscosity  and  conductivity 
coefficients : 


u '  w ' 


-K 


In 

M*z 


T'"'  -  -KmE  <"> 

^  -  -%H 


Above  the  mixed  layer,  turbulent  mixing  makes  such  a  method 
which  assumes  and  constant  unrealistic  and  inaccurate. 
Thus  the  Garwccd  (1977)  turbulence  closure  model  is  invoked 
using  bulk  turbulent  kinetic  energy  equations. 
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The  entrainment  buoyancy  flux  is  given  by: 

$)/h  (12) 

where  <w*  w  •>  and  <E>  are  the  vertical  and  total  components 
respectively  cf  the  turbulent  kinetic  energy.  The  mixed 
layer  total  turbulent  kinetic  energy  equation  and  the  verti¬ 
cal  component  of  turbulent  kinetic  energy  equation  are: 

i_jkiOJ=  ^  -  *gh(^)^/(2Ri*)  -  (<S>*  +  fh)<2>  03) 


<kgh((w’T’).h  “  (w»T')0)/2 

+  (<e>  -  3<^r,))<2>i 

-  (<E>*  ♦  fh)<E>/3 


(H) 


where  u*  =  ( 


is  the  friction  velocity  and 


Ri*=  _ «*  ghflg 

dU  +dV 


is  the  bulk  Richardson  Number. 


The  downward 
puted  as: 


fluxes  of  momentum  -  (CrnF7 )  and  -(v'w‘) 
-(uTwT)_h  =  WgAu  =  - (w'T* ) _h(^j 

-(v'w’)_h  =  w  gAv  =  -(WT’).^ 


are  com- 


(15) 
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where  wg  is  the  entrainment  velocity: 


w 


e 


-(w«T’)_h 

AT 


(16) 


the  entrainment,  heat  flux  (w’T1)  is  the  lower  boundary 

-ri 

condition  for  the  mixed  layer  temperature  profile,  and  a  new 
mixed  layer  depth  and  new  momentum  and  temperature  profiles 
are  calculated.  The  mixed  layer  depth  is  independent  cf  the 
model  level  depths.  If  the  mixed  layer  shallows  rather  than 
deepens  (as  wculd  occur  if  there  were  insufficient  vertical 
turbulent  energy  to  transport  heat  down  to  the  depth  of  the 
existing  mixed  layer),  then  the  depth  where  (w'T* ) — »  0 
becomes  the  new  mixed  layer  depth.  During  mixed  layer  shal¬ 
lowing,  heat,  momentum  and  potential  energy  are  conserved. 
To  prevent  the  levels  _eneath  the  mixed  layer  from  becoming 
unstable  in  the  numerical  process,  a  dynamic  stability  con¬ 
dition  is  imposed  such  that  the  Gradient  Richardson  Number 
is  always  greater  than  1/4. 

The  embedding  of  the  mixed  layer  model  within  the  gen¬ 
eral  circulation  model  requires  communication  between  the 
two.  The  general  circulation  model  is  the  dynamic  portion 
which  calculates  changes  of  u,  v,  and  T  at  each  depth  level 
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due  to  advective  and  diffusive  processes.  These  values  ar; 
then  acted  upcn  fcy  the  mixed  layer  model,  which  in  turn  cal¬ 
culates  the  changes  due  tc  surface  fluxes  and  entrainment 
mixing,  and  calculates  the  new  mixed  layer  depth.  These 
values  are  then  used  ty  the  general  circulation  model. 


III.  SOD  EL  RE  SOLIS  FOR  FRONT  ±  CASE  I 
A.  ANALYSIS  AND  DISCUSSION 

Although  the  front  is  initially  in  geostroFhic  balance, 
the  imposition  cf  atmospheric  forcing  causes  the  front  to 
undergo  a  variety  of  adjustments.  In  Case  I,  a  positive 
wind  stress  cf  1  dyna/cm^  is  applied.  The  front  and  the 
u-comFonent  field  undergo  diffusion,  thereby  giving  rise  to 
an  ageostrophic  u-component,  which  in  turn  creates  a  v-ccm- 
ponent.  As  diffusion  occurs  through  mixing,  and  as  the  sur¬ 
face  expression  of  the  front  is  advected  to  the  left  in 
response  to  the  wind  stress,  the  along-front  jet  decreases 
in  magnitude  as  well  as  in  horizontal  extent.  A  small 
countercurrent  is  present  at  hour  24  due  to  a  reversal  in 
the  slope  of  the  isotherms  in  the  vicinity  of  the  front. 
This  effect  disappears  by  hour  36,  but  then  reappears  at 
hour  48. 

The  base  cf  the  front  on  the  right  side  deepens  (as  evi¬ 
dent  by  the  isotherms  becoming  vertical  to  greater  depths) 
from  about  30  s  initially  to  about  50  m  as  a  result  of  the 
effect  of  turbulent  mixing  and  the  creation  of  ageostrophic 
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velocities  while  undergoing  dynamic  adjustment  (see  Fig.  7). 
An  ageostrcphic  u-component  is  evident  in  Fig.  8  in  the  area 
of  y*53  km  to  61  km.  The  mixed  layer  depth  increases  as  the 
upward  sloping  thermccline  erodes  within  the  front  because 
water  is  transported  across  the  front,  right  to  left.  The 
mixing  of  the  coder  water  with  the  warmer  water  causes  the 
mixed  layer  depth  to  deepen  on  the  left  (warm)  side.  Tha 
axis  of  maximum  v- values  lies  just  above  the  bottom  of  the 
mixed  layer  near  y=5 3  km.  If  the  v-profile  is  examined 
(Fig.  9)  ,  it  is  seen  that  dv/dy  >  0  on  the  right  (cold)  side 
of  the  frontal  axis,  and  dv/dy  <  0  on  tha  left  side.  The 
gradients  cf  v  are  much  stronger  on  the  left  side  of  the 
front  than  they  are  on  the  right.  The  maximum  v-velocity  is 
-8  tc  -10  cm/c.  From  continuity,  dw/dz  <  0  (downwelling)  on 
the  right  side  of  the  axis  and  dw,  )z  >  0  (upwelling)  in  a 
thin  layer  on  the  left  side  that  closely  follows  the  mixed 
layer  depth.  It  is  noted  also  that  the  mixed  layer  depth 
has  deepened  at  the  right  side  of  the  front.  This  phenom¬ 
enon  is  present  throughout  the  model  integration.  Though 
such  a  conditict  was  not  anticipated  prior  to  the  running  of 
the  model  (we  expected  an  upwelling  in  the  mixed  layer  on 
the  right  side  as  a  result  of  the  mass  transport  over  the 
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front  towards  the  left) ,  this  depression  in  depth  is 
believed  tc  be  physically  realistic.  It  is  suspected  that 
this  increase  in  the  nixed  layer  depth  occurs  in  response  to 
no n-staticnary  adjustment  to  the  atmospheric  forcing,  as  it 
also  occurs  in  Case  II  where  the  oppositely- directed  wind 
stress  is  applied.  Because  the  mixed  layer  depth  responds 
to  the  vertical  integral  of  transport,  the  fact  that  the 
mixed  layer  is  shallower  within  the  front  than  in  the  far 
field  blocks  the  flow  frcm  penetrating  to  depth.  As  the 
water  encounters  this  "obstacle",  it  "backs  up"  on  the  right 
side  and  thus  fcrcas  a  depression  in  the  mixed  layer.  This 
explanation  must  be  tempered  by  the  fact  that  these  figures 
are  instantaneous  snapshots  in  time,  and  can  be  deceiving  in 
attempting  tc  describe  processes  that  have  been  evolving 
over  a  length  of  time.  Another  possibility  is  that  soma 
type  of  internal  standing  wave  has  been  excited  which  is 
causing  fluctuations  in  the  mixed  layer  depth. 

The  influence  of  the  ageostrophic  u-velocity,  the  core 
of  the  v-velocity  lying  along  the  well-mixed  layer  a  the 
depression  of  the  mixed  layer  on  the  right  side  of  the  front 
is  felt  at  depths  greater  than  100  m.  Upwelling  and  down- 
welling  are  evident  at  these  depths  (Fig.  7)  beneath  the 
right  edge  of  the  front. 
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Figure  8.  Front  1  Case  I  12-Hour  U-Velocity.  Contour 
interval  is  8  ca/s . 
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Interesting  patterns  in  the  v  profiles  evolve 


at  hours 


24  and  48.  At  24  hours  a  small  positive  v  exists.  This  is 
contrary  to  the  expected  Ekman  transport  direction  in  the 
surface  layer,  and  is  probably  the  result  of  an  oscillation 
created  by  the  inertial  action.  Furthermore,  at  both  24  and 
48  hours,  two  cores  of  a  positive  v-component  are  present 
symmetrically  centered  about  the  base  of  the  mixed  layer  on 
either  side  cf  the  front.  These  are  not  present  in  the 
front,  but  fora  outside  the  frontal  area  and  extend  to  the 
far  field  where  the  mixed  layer  is  relatively  stable. 

The  across-front  velocity  component  is  in  consonance 
with  the  direction  of  the  net  Ekman  transport  at  all  times 
except  hour  24.  As  noted,  this  velocity  is  at  a  maximum  at 
the  base  of  the  mixed  layer  rather  than  at  the  surface.  It 
is  ncn-zerc  beneath  the  mixed  layer  in  the  middle  of  the 
front  as  well.  Thus,  it  appears  that  interfacial  mixing 
does  occur  in  this  case,  and  it  may  be  contributing  to  the 
vertical  mixing  cf  heat  in  the  thermocline  at  this  location. 

B.  SQMHAfll  —  CASE  I 

Under  the  influence  of  positive  wind  stress  and  no  heat¬ 
ing,  the  frcct  has  diffused  and  has  been  advected  in  the 
direction  of  the  surface  Ekman  transport.  Because  the  frcnt 


is  unsteady  under  an  applied  wind  stress,  ageostrcphic 
velocities  are  created.  Fart  of  the  response  is  frictional 
due  to  the  turbulent  boundary  layer  or  aixad  layer  pro¬ 
cesses,  and  part  is  inertial  due  to  the  applied  wind  stress 
boundary  condition.  &  noticeable  deepening  of  the  mixed 
layer  occurs  at  the  front  on  the  right  (cold)  side,  and  the 
horizontal  variability  in  nixed  layer  depth  and  temperature 
are  reduced  significantly  under  the  influence  of  both  verti¬ 
cal  and  cross-frontal  mixing. 
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IV.  KCDEL  SESOLTS  FOE  FRONT  J.  CASE  II 
A.  ANALYSIS  AND  CISCOSSICN 

In  this  case,  a  negative  wind  stress  (directed  into  the 
page)  is  applied.  Wind  stirring  oust  overcome  static  sta¬ 
bility  as  the  warmer  water  will  override  the  cooler  water. 
It  will  be  seen  that  the  wind  stress  of  1  dyne/ca^  is  suffi¬ 
cient  to  mix  the  warmer  water  which  is  transported  toward 
the  ccoler  water.  There  is  no  evidence  of  stratification  of 
shallow  waters.  Recall  from  Case  I  that,  as  the  denser 
water  from  the  right  side  was  being  transported  to  the  left, 
it  had  to  overcome  static  stability  in  order  to  upwell  and 
override  the  front,  and  that  the  nixed  layer  depth  profile 
posed  a  vertical  barrier  tc  the  transport.  In  this  casa, 
the  base  cf  the  mixed  layer  is  at  an  incline  in  the  direc¬ 
tion  of  transport  and  static  stability  need  not  be  overcome 
in  the  transport  process.  This  leads  to  the  expectation 
that  cross-frontal  mixing  may  be  minimal  in  this  case. 

Diffusion  in  the  temperature  field  again  is  evident  over 
the  48  hours,  but  this  ti«e  it  occurs  primarily  at  the  edges 
of  the  front.  A  strong  horizontal  temperature  gradient  of 
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0.6  deg  C/1.2  km  in  the  aiddle  of  the  front  is  maintained  up 
through  48  hours  (see  Pig.  10  as  representative  of  this) 
which  was  not  present  in  Case  I.  This  lends  credence  to  the 
hypothesis  that  cross-frontal  mixing  is  minimal. 

Noticeable  mixing  occurs  on  the  cold  side  of  the  frcnt. 
The  isotherms  are  vertical  to  greater  depths  than  they  were 
in  Case  I  (see  Fig.  10) .  The  turbulent  mixing  alone  net 
only  overcomes  static  stability  but  is  also  able  to  mix  to  a 
relatively  greater  depth  cn  the  right  side  of  the  front. 
The  water  that  is  being  transported  across  the  front  is 
apparently  immediately  mixed  by  wind  stirring  as  no  plume  of 
the  warm,  less  dense  water  is  present  cn  the  right  side. 

The  alcng-front  current  again  forms  an  ageostrophic  com¬ 
ponent  as  it  undergoes  adjustment,  and  is  especially  evident 
on  the  right  side  of  the  front  around  y=55  km  (Fig.  11) 
where  it  exists  to  a  greater  extent  both  horizontally  and 
vertically  than  was  seen  in  Case  I.  The  core  of  the  alcng- 
front  jet  dees  decrease  in  speed  and  diffuse  in  time.  By 
hour  36  the  maximum  value  of  u  has  been  reduced  to  slightly 
more  than  8  cm/s.  The  creation  of  this  ageostrophic  veloc¬ 
ity  accompanied  by  a  v-compcnent  leads  to  strong  dovnwelling 
around  y=55  km  underneath  the  front  (Fig.  10)  and  remains 
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over  the  48  hour  integration.  A  strong  core  of  negative 
across-front  velocity  of  over  -6  cm/s  is  created  in  this 
area  (Fig.  12),  and  there  is  an  attendant  sharp  upwelling 
in  the  isotherms  at  75  m  and  vave-like  fluctuations  present 
to  depths  beycnd  175  m.  In  the  areas  where  dv/dy  >  0  (dw/dz 
<  0) ,  downwelling  is  seen  in  the  front,  and  where  dv/dy  <  0 
(dw/dz  >  0) ,  upwelling  is  seen.  In  the  surface  layer,  the 
v-component  is  positive  and  in  the  direction  of  Ekman  trans¬ 
port.  It  is  concentrated  in  a  jet  in  the  area  of  the  tight¬ 
est  horizontal  temperature  gradient  over  the  shallow  peak  in 
the  mixed  layer  depth.  All  of  these  tendencies  and  trends 
persist  over  the  48  hour  integration. 

A  deepening  in  the  mixed  layer  occurs  to  the  right  of 
the  front  near  y=55  km  and  progresses  outward  in  time.  The 
shallow  peak  in  the  mixed  layer  depth  within  the  front  is 
mostly  preserved,  though  some  erosion  and  mixing  obviously 
has  occurred  (compare  Figs.  13  and  14).  This  further  sup¬ 
ports  the  hypothesis  that  minimal  cross-frontal  mixing 
occurs.  Examination  of  the  corresponding  v-fields  also  show 
little  cross-frontal  nixing. 
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B.  SOMMABY  —  CASE  II 

In  this  case  of  a  forcing  of  a  negative  wind  stress  and 
no  heating,  the  temperature  structure  of  the  front  retains  a 
strong  horizontal  gradient  over  the  48-hour  integration. 
Diffusion  cf  the  front  occurs,  but  it  is  concentrated  at  the 
boundaries  of  the  front.  Sind  stirring  overcomes  the  static 
stability  cf  transported  water  once  it  crosses  the  frcnt. 
The  creaticr  under  adjustment  of  a  large  ageostrophic  u-ve- 
locity  and  an  accompanying  v-velocity  give  rise  to  strong 
up  welling  and  dcwnwelling  features  within  the  front.  The 
shallow  peak  cf  the  mixed  layer  depth  within  the  front  is 
preserved  ever  time.  The  fact  that  this  peak  is  maintained, 
the  horizontal  temperature  gradient  retains  a  tight  struc¬ 
ture,  and  the  v-field  is  mainly  constrained  to  the  agecs- 
trophic  areas  and  surface  layer,  suggests  that  intarfacial 
or  cress- f  rental  mixing  is  minimal. 
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Figure  11.  Frcnt  1  Case  II  36-Hour  0-Velocity.  Contour 
interval  is  4  cn/s. 
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Figure  12.  Front  1  Case  II  12-Hour  V-Velocity.  Contour 
interval  2  cb/s. 
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H  at  hour  48 


V.  aCDEL  RE  SOUS  FOR  FRONT  1  CASE  III 


A.  ANALYSIS  AND  DISCUSSION 

A  uniform  surface  heating  is  applied  with  the  identical 
wind  stress  cf  Case  I.  The  magnitude  of  each  forcing  is 
such  that  influences  from  both  are  expected  to  be  mani¬ 
fested.  In  examining  Figs.  15  and  16,  the  obvious  effect  of 
heating  is  the  adjustment  of  the  temperature  structure  on 
either  side  cf  the  front.  At  first  examination,  it  appears 
that  the  heating  is  having  the  opposite  effect  of  what  is 
expected  -  that  is,  it  appears  to  cool  the  mixed  layer 
rather  than  warm  it.  In  reality,  the  heat  is  carried  dcvn 
to  depths  cf  50  m  and  warms  the  water  there.  Cooling  occurs 
between  32  and  50  m  depth,  which  is  readily  seen  in  Fig.  17. 
Hiring  may  create  a  dynamic  instability  in  the  thermocline, 
and  vertical  diffusion  of  heat  by  Ryed^T/dy^  may  cause  the 
intermediate  water  to  be  ceded  and  the  deep  mixed  layer  za 
be  warmed.  This  basic  structure  remains  over  the  48-hcur 
integration.  The  frent  undergoes  diffusion  and  is  advected 
to  the  left,  as  occurred  in  Case  I.  The  front  is  manifested 
to  a  lesser  depth  in  accordance  with  the  heating  process  - 
downwelling  and  mixing  have  been  moderated  by  the  heating. 
Though  some  wave-like  motion  is  seen  in  the  temperature 
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field  beneath  the  mixed  layer,  the  magnitude  is  not  as  great 
as  in  the  case  of  wind  stress  alone.  The  mixed  layer  depth 
has  shallowed  to  about  43  m  under  heating  alone  in  the  areas 
outside  the  front  (see  Fig.  16).  This  shallowing  occurs 
because  the  atacspheric  forcing  (heating  and  wind)  is  input 
as  a  body  force  ever  the  entire  mixed  layer.  The  heating 
has  the  effect  of  making  the  mixed  layer  depth  more  uniform. 
The  marked  deepening  in  the  mixed  layer  depth  that  occurred 
in  Cases  I  and  II  at  the  right  (cold)  side  of  the  front  is 
not  present  in  this  case.  The  deepening  of  the  mixed  layer 
affects  the  u-field.  The  ageostrophic  u-component  does  not 
develop  to  as  great  an  extent  as  it  did  ir.  Casa  I.  The  gen* 
eral  along-frent  current  appears  much  as  it  did  in  Case  I, 
except  that  its  vertical  extent  is  also  constrained  by  the 
shallower  mixed  layer  depth.  In  both  cases,  maximum  speeds 
agree  over  the  48-hour  integration,  and  they  are  advacted 
with  the  front  as  it  responds  to  the  mass  transport. 

The  v-field  follows  the  same  trends  as  it  did  under  wind 
stress  alone.  At  12  hours,  v  is  directed  to  the  left  as 
expected  following  the  transport,  and  at  24  hours  the  v-ve- 
locity  reverses  and  establishes  the  layers  of  concentrated 
positive  v-values  centered  symmetrically  about  the  mixed 
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Figure  16.  Front  1  Case  III  12- Hour  Mixed  Layer  Depth. 
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Figure  17.  Front  1  Case  III  -  T  vs  2  Profile.  Profile 
taken  at  y»35  ka  for  t=0  and  t*12  hours. 


layer  depth.  At  36  hairs,  v  returns  to  a  negative  value  ana 
at  48  hours,  though  still  negative,  v  is  such  decreased  in 
magnitude.  Again,  cells  cf  positive  v-velocities  are  pres¬ 
ent  at  about  the  mixed  layer  depth.  The  across-front  veloc¬ 
ity  does  penetrate  the  front  and  the  shallow  peak  of  the 
mixed  layer  (see  Fig.  18)  which  is  indicative  of  cross-frcn- 
tal  mixing.  The  shallow  section  of  mixed  layer  itself  is 
eroded  and  smoothes  cut  ever  the  48-hour  integration  as  it 
did  in  Case  I  (Fig.  19). 

B.  SUMHAEY  --  CASE  III 

The  effects  cf  both  wind  stress  and  heating  are  clear  in 
this  case.  The  wind  stress  induces  a  mass  transport  to  the 
right,  moving  cooler,  denser  water  into  warmer,  less  dense 
water.  This  water  not  only  apwells  and  moves  over  the 
front,  but  also  moves  across  the  frontal  interface.  At  24 
and  48  hours,  the  v- field  appears  to  be  influenced  by 
inertial  oscillations  with  a  reversal  in  direction  occurring 
in  the  surface  layer  at  24  hours  and  a  near-reversal  at  48 
hours.  Cores  cf  positive  velocities  are  established  cen¬ 
tered  about  the  mixed  layer  depth  outside  of  the  front.  The 
front  both  diffuses  and  moves  to  the  left  in  time. 
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Pigurs  19.  Prcnt  1  Case  III  48- Hour  Mixed  Layer  Depth. 
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Heating  shallows  the  nixed  layer  depth  as  expected,  and 
thereby  restricts  the  vertical  extent  of  the  velocity 
fields.  Heating  also  prevents  the  immediate  deepening  of 
the  fixed  layer  on  the  right  side  of  the  front,  which  was 
also  seen  in  both  Cases  I  and  II.  The  heat  manifests  itself 
at  the  initial  mixed  layer  depth  where  the  isotherms  become 
more  stratified. 
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VI.  MODEL  RESULTS  FO  R  FROKT  1  CASE  IV 

A.  ANALYSIS  AND  DISCUSSION 

Case  IV  combines  a  negative  wind  stress  and  a  uniform 
surface  heating.  The  mixed  layer  shallows  noticeably  and 
heat  is  carried  to  the  bottom  of  the  mixed  layer  at  depths 
of  55  xo  70  a  as  in  Case  III  (see  Figs.  20  and  21).  The 
heating  dees  have  an  interesting  effect  when  compared  with 
the  wind- strass-cnly  forcing  of  Case  II.  The  sharp  dip  ir. 
the  isotheras  in  the  middle  of  the  front  is  not  evident  here 
as  it  was  in  Case  II.  Horizontal  diffusion  of  the  tempera¬ 
ture  structure  occurs  mors  rapidly  and  attains  a  gradient  of 
1.6  deg  C/5.6  km  by  hour  12  (see  Fig.  20)  .  The  along-fronx 
velocity  is  limited  in  its  vertical  development  on  the  right 
side  of  xhe  front  near  y=57  km  (Fig.  22)  and  follows  both 
the  horizontal  and  verxical  diffusion  tendencies  of  the  den¬ 
sity  front.  The  across-front  velocity  has  a  narrow,  verti¬ 
cally-oriented  core  of  negative  values  with  a  maximum  of  -6 
cm/s  cn  the  right  side  of  the  front  between  25  and  50  m 
depth  as  in  Case  II.  Since  this  maximum  of  v-velocity 
exists  in  an  area  where  there  is  little  u-velocity,  its 
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existence  cannot  be  attributed  to  a  response  to  the  agscs- 
trophic  u-ccmpcnent.  The  position  of  this  v-velocity  corre¬ 
lates  well  with  the  upwelling/downwelling  pattern  evident  in 
the  isotherms  in  Pig.  20,  and  also  with  wave  patterns  in  the 
isotherms  over  the  entire  depth  extent  in  the  right  half  of 
p'iq.  23.  We  suspect  that  such  a  large  area  of  v-coaponent 
created  below  the  nixed  layer  is  a  response  to  the  warming 
of  the  water  at  the  initial  nixed  layer  depth  which  has 
shoaled,  re-oriented  the  isot haras  and  created  vertical 
velocity .  Furthermore,  wind-generated  mixing  forces  verti¬ 
cal  motions  in  the  water  nass  to  a  great  dapth,  inducing  a 
v-conponent  there. 

In  this  case,  the  front  diffuses  horizontally  more  rap¬ 
idly  than  under  wind  stress  forcing  alone.  The  front 
spreads  like  a  warn,  less-der.se  plume  overriding  the  cool, 
denser  water  (compare  Fig.  24  to  Fig.  20)  .  The  along-frcnt 
current  shifts  with  the  front.  The  mixed  layer  loses  its 
shallow  section  within  the  front  and  smoothes  out,  which  did 
not  happen  in  Case  II.  The  peak  in  the  mixed  layer  was 
maintained  much  longer.  With  heating,  there  is  less  of  a 
temperature  difference  and  less  of  a  shallowing  from  the 
far-field  mixed  layers  to  the  shallow  mid-frontal  mixed 
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Figure  20.  Front  1  Case  IV  12-Hcur  Tenperature .  Contour 
interval  is  0.2  deg  C. 


68 


U  at  hour  12 


K 


Depth  (m) 


V  at  hour  12 


Figure  23.  Front  1  Case  IV  12-Hcur  V-Velociry.  contour 
interval  is  2  ca/s. 
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layer,  which  probably  accounts  for  the  quicker  breakdown  of 
the  shallow  section  of  the  mixed  layer. 

The  v-field  in  the  surface  layer  is  advectsa  to  the 
right,  with  a  maximum  value  located  at  the  surface  within 
the  front.  By  hour  24,  strong  horizontal  cells  of  a  v-ccm- 
ponect  have  formed  about  the  mixed  layer  depth.  In  this 
case,  the  values  are  negative,  whereas  in  Cases  I  and  III 
they  were  positive.  (Case  II  never  formed  them.).  These 
cells  disappear  at  hour  36  and  reappear  at  hour  48,  as  they 
did  in  Cases  I  and  III.  The  mechanism  that  explains  the 
appearance  of  negative  v  values  is  unclear. 

The  trends  to  48  hours  show  further  diffusion  and  advec- 
tion  of  the  warm  frontal  plume  to  the  right.  There  is 
smoothing  of  the  horizontal  distribution  of  mixed  layer 
depth.  The  surface  front  patterns  of  u-  and  v-velocity  are 

all  aavected  to  the  right. 
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B.  SOMMABY  —  CJ>SS  IV 

Application  cf  a  negative  wind  stress  and  surface  heat¬ 
ing  results  in  the  front  moving  as  a  warm  overriding  plume 
to  the  right.  The  mixed  layer  initially  shsllcws  due  to  the 
heating,  and  ever  the  48  hour  integration  breaks  down  and 
smoothes  cut  in  the  horizontal.  The  along-front  velocity  is 
advected  and  diffused  in  the  mixed  layer.  The  pattern  cf 
across-frcnt  velocity  shifts  to  the  right  in  the  surface 
layer,  with  its  maximum  value  centered  over  the  shallowest 
portion  of  the  mixed  layer.  A  v-ccmponant  develops  to  quite 
a  great  depth  on  the  right  side  of  front  due  to  motion  gen¬ 
erated  by  heating  and  turbulent  mixing.  Cells  of  negative 
v-component  are  formed  at  hours  24  and  48  which  are  centered 
about  the  mixed  layer  depth. 
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VII.  MODEL  RESULTS  FOR  FRONT  2  CASE  I 

A.  ANALYSIS  AND  DISCOSSICN 

The  initial  conditions  for  Front  2  are  shown  in  Figs.  4, 
5,  and  6.  The  front  extends  throughout  the  depth  of  the 
basin  and  is  oriented  such  that  warmer,  less  dense  water  is 
on  the  right  and  cooler,  denser  water  is  on  the  left.  Many 
of  the  dynamical  features  of  Pront  1  are  expected  to  be 
repeated  fcr  Front  2.  In  comparing  Front  2  to  Front  1,  the 
horizontal  axes  are  oppositely-directed  to  the  frontal  ori¬ 
entation.  The  thermocline  slopes  upward  in  the  positive 
y-directicn  for  Front  1,  tut  downward  for  Front  2. 

The  wind  stress  cf  Case  I  applied  to  Front  2  creates  a 
transport  directed  to  the  left.  Warmer,  less  dense  water 
moves  towards  the  cooler,  denser  water.  In  this  case,  dif¬ 
fusion  in  the  temperature  field  is  present  for  the  duration 
of  the  integration.  The  horizontal  temperature  gradient 
remains  large  (1.C  deg  C/1.2  km)  in  the  middle  of  the  front 
with  the  majority  of  the  diffusion  occurring  at  the  right 
boundary.  Mixing  on  the  cool  (left)  side  of  the  front 
extends  from  the  surface  to  1 35  a  which  is  below  tfct-  initial 
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depth  of  the  nixed  layer.  Results  of  downwelling  are  seen 
at  the  left  frontal  boundary,  near  y=42  ha.  Apparently, 
there  is  no  chance  for  a  stable,  overriding  warm  plume  to 
develop  before  it  is  nixed  into  the  ambient  cool  water  peel. 
These  phenomena  in  the  temperature  field  were  also  seen  in 
Case  II  for  Front  1.  Much  more  dramatic  oscillations  occur 
in  the  isotherms,  primarily  on  the  right  side  between  53  km 
and  70  km,  where  the  entire  field  is  inclined.  They  are 
also  visible  cr  the  left  side  superimposed  upon  an  inclina¬ 
tion  of  the  entire  field  sloping  downward  to  the  right  (Fig. 
25).  The  isctherms  show  a  definite  trend  towards  becoming 
more  uniformly  itclined  beneath  the  mixed  layer  across  the 
entire  horizontal  expanse,  with  the  steeply-inclined  field 
on  the  right  becoming  more  horizontal.  Downwelling  causes 
the  left  side  to  slope  downward  to  the  right  and  upwelling 
lifts  the  isotherms  on  the  right  side. 

The  mixed  layer  depth  undergoes  major  changes  over  the 
entire  domain.  Cn  the  left  side,  it  deepens  to  where  the 
depth  in  the  far  field  is  approximately  125  m  (Fig.  26)  . 
Since  this  deepening  occurs  by  3  hours  of  integration,  such 
a  drastic  deepening  probably  resulted  from  the  initial  mixed 
layer  depth  being  inserted  at  too  shallow  a  depth  in  this 
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T  at  hour  36 


Figure  25.  Frcnt  2  Case  I  36-Hour  Tenperature. 
interval  is  0.5  deg  C. 


Contour 


area.  As  in  Front  1  Case  II,  a  downward  bulge  in  the  mixed 
layer  erterdirg  to  155  m  occurs  at  the  left  frontal  boundary 
at  y*43  km  (Fig.  26)  .  In  this  case,  the  downward  bulge  does 
not  correspond  to  rhe  side  on  which  the  waters  would  tend  to 
"back  up"  prior  to  crossing  over  the  front  as  in  Front  1 
Case  I.  As  strong  nixing  is  implied  by  the  vertical  orien¬ 
tation  of  the  isotherms  at  that  location,  the  bulge  is 
believed  to  te  associated  with  a  cooling  of  the  waters  in 
that  area  as  a  result  of  enhanced  vertical  nixing  caused  by 
the  ageostrophic  along-front  velocity  shear  (Fig.  27) . 
Within  the  front,  a  shallow  peak  in  the  mixed  layer  depth  is 
maintained  at  the  left  boundary,  though  over  time  it  does 
deepen  slightly  and  erode  through  mixing.  There  is  much 
variability  ir.  the  mixed  layer  depths  on  the  right  side  of 
the  front,  but  it  appears  that  the  general  trend  is  for 
deepening  to  occur  in  the  far  field  beyond  y=55  km  (Fig. 
26)  . 

Along  with  the  isotherms,  the  along-front  u-velocity 
diffuses  and  decreases  its  maximum  speed  over  time.  As 
already  seen,  an  ageostrophic  component  forms  on  the  left 
frontal  boundary  and  stimulates  mixing  there.  The  remnants 
of  a  tight  horizontal  velocity  gradient  are  seen  after  48 
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H  at  hour  36 


Figure  26.  Front  2  Case  I  36-Hour  Hired  Layer  Depth 


U  at  hour  12 


Y-distance  (km) 


Figure  27.  Front  2  Case  I  12-Hour  0-Velocity.  Contour 
"  interval  is  1C  cm/s. 


hours,  though  little  cr  nc  advecticn  of  the  along-front  cur¬ 
rent  occurs  in  response  to  the  wind-generated  transport. 
Have-like  oscillations  appear  to  propagate  along  the  deeper 
extremities  cf  the  u-field,  but  no  phase  relation  between 
isotachs  cr  with  ether  fields  is  readily  apparent. 

The  acrcss-front  v-velocity  does  correlate  with  the 
direction  cf  transport  up  through  hour  9.  An  oscillatory 
pattern  exists  iE  the  v  field  over  the  48-hour  integration. 
At  hour  12,  the  entire  v-field  is  positive,  but  this  may  be 
a  response  to  the  agecstrcphic  u-velocity,  and  it  is  strong 
enough  to  overcome  the  acrcss-frcnt  transport  velocity.  At 
24  hours,  v  is  in  the  direction  cf  transport  and  a  maximum 
of  -20  cm/s  appears  at  the  surface  within  the  front  over  the 
shallowest  section  of  mixed  layer  depths  at  y=43  km  (Fig. 
28) ,  and  remains  through  hour  36 .  This  is  reassuring  as  it 
was  expected  that  the  inherent  static  stability  of  the  water 
masses  would  uiniaize  deep  cross-frontal  mixing,  similar  to 
what  was  seen  in  Front  1  Case  II.  The  isotachs  very  defi¬ 
nitely  terminate  at  the  left  boundary  of  the  front,  but  They 
also  exist  beneath  the  mixed  layer  and  do  cross  the  front 
beneath  the  mixed  layer.  At  hour  48,  v  is  positive  again, 
perhaps  due  to  inertial  motions.  It  is  proposed  (though  it 
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is  certainly  net  as  clear  as  it  was  for  Front  1)  that 
cr oss-frontal  nixing,  though  inhibited,  occurs  here  to  a 
greater  extent  than  it  did  for  Pront  1.  The  majority  of  the 
mass  transport  appears  to  exist  in  the  mixed  layer  and 
passes  over,  rather  than  through,  the  shallow  mid-frontal 
depth  of  the  nixed  layer.  The  reasons  for  the  variations 
between  the  fronts  may  be  due  to  the  greater  u-component 
(both  in  speed  and  in  areal  coverage)  present  with  Front  2. 

B.  SUHMABI  --  CASE  I 

Several  complexities  have  been  incorporated  into  the 
modeling  of  Front  2.  First,  the  isotherms  slope  downward  to 
the  right  ever  the  entire  tight  half  of  the  basin.  Hence, 
the  front  not  cnly  exists  at  the  surface  but  is  "connected" 
through  the  depth  of  the  tasin.  This  establishes  an  initial 
u-velocity  which  exists  over  the  same  area.  Adjustment  to 
the  wind  stress  may  not  be  confined  solely  to  the  upper  few 
meters,  as  dynaiic  responses  may  be  transmitted  to  depths 
beyond  those  directly  affected  in  the  surface  layer.  Also, 
the  mixed  layer  depth  is  artificially  inserted  as  an  initial 
condition,  and  nay  net  be  the  most  perfect  fit  to  the  temp¬ 
erature  field.  Though  this  mixed  layer  configuration  may 
exaggerate  cr  suppress  actual  features,  it  is  nevertheless 
felt  to  be  a  reasonable  approximation  of  an  actual  front. 
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The  responses  are  much  more  complicated  in  structure  in 


Front  2  than  they  were  in  Front  1.  Beyond  12  hours,  the 
front  shows  little  advective  response  to  the  mass  transport. 
One  reason  this  occurs  is  because  the  front  exists  over  such 
a  great  depth.  The  original  thermal  structure  and  tempera¬ 
ture  gradient  are  maintained  in  the  middle  of  the  front, 
with  horizcntal  diffusion  occurring  primarily  at  the  right 
boundary.  Mixing  in  the  surface  layer  occurs  as  the  isot¬ 
herms  become  vertical  to  greater  depths  over  time.  At  the 
laft  boundary  cf  the  front,  downwelling  is  quite  evident  and 
the  waters  mix  easily  due  tc  the  creation  of  an  ageostrophic 
u-component  which  enhances  the  wind  stirring.  Wave-like 
oscillations  are  evident  in  the  temperature  field  beneath 
the  mixed  layer,  and  the  entire  deep  temperature  field 
becomes  mere  uniform  in  slope  over  the  48-hour  integration. 

Tha  mixed  layer  depth  must  be  considered  in  three  sepa¬ 
rate  regimes  -  left  side,  frontal,  and  right  side.  On  the 
left,  the  mixed  layer  depth  increases  by  almost  50  m  at  the 
frontal  boundary,  where  mixing  is  enhanced  due  to  the  influ¬ 
ence  of  the  acecstrophic  current.  The  mixed  layer  bulges 
downward,  then  rises  sharply  to  its  shallow  peak  at  the 
front.  This  shallow  section  remains  throughout  the  48  hour 
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integration  though  it  dees  deepen  slightly  due  to  mixing. 
On  the  right  side,  the  mixed  layer  deepens  and  becomes  more 
horizontal  though  there  are  strong  fluctuations  over  inter¬ 
vals  of  4  km  and  less. 

Two  primary  trends  are  observed  in  the  along-front  u-ve- 
locity:  1)  the  u-component  diffuses  in  time,  still  retain¬ 
ing  its  basic  shape  but  decreasing  in  speed  and  horizontal 
gradient;  2)  an  ageostrcphic  component  is  evident,  prima¬ 
rily  extending  outwards  f  icm  the  left  boundary  of  the  front 
and  which  is  believed  to  '‘drive"  the  downward  bulge  in  the 
mixed  layer  there. 

The  across-front  velocity  component  is  concentrated  at 
the  surface  in  the  direction  of  transport  between  0  and  9 
hours  and  at  24  and  36  hours.  The  front  poses  no  barrier  to 
the  v-velocity  at  12,  24,  or  48  hours.  The  36-hour  fields 
are  the  crly  cr.es  which  have  essentially  no  cross-front 
velocity  component  other  than  in  the  surface  layer.  Because 
the  initial  horizontal  temperature  gradient  is  preserved  in 
the  middle  of  the  front,  and  because  the  shallow  portion  of 
the  maxed  layer  depth  is  still  evident  after  48  hours,  sug¬ 
gest  that  little  cross-f rental  mixing  occurs. 
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fill.  HOD  El  RESULTS  POP  FRONT  2  CASE  II 

A.  ANALYSIS  AND  DISCUSSION 

With  the  wind  stress  opposite  from  that  of  Case  I,  the 
net  Ekman  transport  is  left  to  right,  and  cooler,  denser 
water  will  be  transported  toward  warmer,  less  dense  water. 
There  are  similarities  here  to  the  results  of  Front  1  Case 
I.  The  front  becomes  much  more  diffuse  at  the  surface  and 
the  horizontal  temperature  gradient  is  not  preserved  as  it 
was  in  the  preceding  case  (see  Fig.  29)  .  It  has  decreased 
to  1.5  deg  C/6.4  km.  The  horizontal  field  on  the  left  side 
reorients  by  inclining  down  to  the  right,  and  the  inclined 
field  on  the  right  becomes  more  horizontal.  The  Total 
effecT  is  to  take  the  temperature  field  more  level  beneath 
the  mixed  layer.  Wave-like  fluctuations  are  also  seen  in 
the  isotherms  (Fig.  29) .  Vertical  mixing  can  be  readily 
seen  to  have  occurred  within  the  front  by  the  steep  slope  of 
the  isotherms. 

The  mixed  layer  depth  immediately  adjusts  in  the  left 
far  field  tc  arcund  125  m  as  it  did  in  Case  I,  but  a  much 
greater  downward  bulge  in  the  mixed  layer  forms  at  the  left 
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T  at  hour  36 
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Figure  29.  Front  2  Case  II  36-Hour  Temperature.  Contour 
interval  is  0.5  deg  c. 


boundary  cf  the  front  around  y*4  2  km  and  remains  over  the  43 
hour  integration  (Fig.  30).  This  result  is  consistent  with 
the  interpretation  proposed  for  Front  1  Case  I,  where  the 
shallower  mixed  layer  within  the  front  acted  as  a  barrier  to 
mass  transport,  causing  the  transported  water  to  "back  up". 
The  mixed  layer  peak  within  the  front  deepens  and  erodes. 
The  right  side  far  field  is  again  variable  but  shows  the 
same  tendencies  to  deepen. 

The  agecstrophic  u-velccity  (not  shown)  is  very  similar 
to  the  isctach  pattern  which  evolved  for  Case  I.  A  cere  of 
maximum  velocity  is  advected  no  the  right  which  correlates 
with  the  advecticn  of  the  front.  The  biggest  difference  is 
that  the  isotachs  appear  at  a  greater  depth  just  inside  the 
left  frontal  boundary. 

The  expected  direction  of  the  across- front  velocity  is 
left  to  right,  the  direction  of  the  net  Ekman  transport.  As 
in  Case  I,  a  wave  oscillation  is  evident.  At  12  and  48 
hours  this  is  clearly  so,  with  the  velocity  field  existent 
along  the  entire  section  cf  the  front  from  the  base  of  the 
front  at  around  150  m  up  to  the  surface  (Fig.  31).  The  max¬ 
imum  velocity  occurs  at  a  depth  beneath  the  shallow  mid- 
frontal  section  of  the  mixed  layer,  which  seems  to 
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contribute  to  cross-frontal  mixing,  erosion  of  the  mixed 
layer  depth  within  the  frcnt,  and  weakening  of  the  frontal 
horizontal  temperature  gradient.  At  24  hours,  the  v-velcc- 
ity  field  is  reversed,  and  negative  velocities  predominate 
in  the  frcrt  and  beneath  the  mixed  layer.  Inertial  motions 
probably  account  for  this  reversal. 

B.  StJHMABY  -  CASE  II 

Ihe  temperature  structure  of  the  front  shows  vertical 
mixing  and  horizontal  diffusion  occurring  over  the  48-hour 
integration.  The  horizontal  temperature  gradient  is 
decreased  due  to  the  effect  of  convective  mixing  of  the 
cooler  water  which  is  transported  to  the  warmer  water  by 
wind  stirring.  The  across-front  velocity  patterns  penetrate 
the  mixed  layer  within  the  front  and  establish  cross-frontal 
mixing.  The  entire  temperature  field  in  the  vicinity  of  the 
front  undergoes  an  adjustment  which  inclines  all  the  isot¬ 
herms  downward  tc  the  right  in  a  rather  uniform  pattern. 

The  mixed  layer  undergoes  major  adjustment  over  the 
48-hcur  integration,  but  the  primary  trends  are  a  deepening 
on  bcth  sides  cf  the  front,  a  leveling  off  on  the  right  side 
and  a  breakdown  cf  the  shallow  section  within  the  front  due 
to  mixing  and  diffusion.  By  extrapolating  the  processes 
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Figure  31.  Frcnt  2  Case  II  48-Hour  V-Velocity. 
interval  is  5  cm/s. 


Contour 


occurring  at  48  hours,  it  is  presumed  that  the  mixed  layer 
depth  would  eventually  even  out  with  little  or  r.o  remains  of 
a  frontal  peak.  As  in  other  cases,  a  marked  deepening  in 
the  nixed  layer  occurs  at  the  left  boundary  of  the  front. 
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IX.  MODEL  REVOLTS  FOR  FRONT  2  CASE  III 
A.  ANALYSIS  AND  DISCUSSION 

Constant  surface  heating  is  combined  with  a  positive 
wind  stress.  As  the  net  Ikman  transport  will  be  from  right 
to  left,  and  with  heating,  we  expect  to  see  the  formation  of 
a  warm  plume  cf  water  on  the  surface  in  the  cool  (left) 
ambient  pool.  The  surface  heating  should  enhance  the  inher¬ 
ent  static  stability  and  prevent  the  wind  from  completely 
mixing  the  transported  warm  waters.  This  effect  does  happen 
prior  to  the  12  hour  point  and  becomes  more  pronounced  over 
the  entire  integration  (see  Fig.  32).  Downwelling  and  mix¬ 
ing  cf  the  warmed  surface  waters  with  the  cooler  underlying 
waters  are  observed  at  the  leading  edge  of  the  progressing 
front.  Upwelling  on  the  right  side  of  the  cooler  water 
which  replaces  water  lost  due  to  transport  is  clearly  seen 
by  comparing  successive  temperature  fields.  The  front  shews 
a  tendency  tc  almost  split  at  30  m  depth  into  a  shallow  sur¬ 
face  frost  and  a  deeper  frent.  Perhaps  a  different  combina¬ 
tion  of  forcing  would  make  this  happen. 
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Figure  32.  Front  2  Case  III  36-Hour  Temperature.  contour 
interval  is  0.2  deg  c. 
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The  isotherms  are  much  smoother  than  for  the  oases  of 
wind  stress  alone.  as  in  Case  I,  the  originally- horizontal 
field  on  the  left  side  becomas  inclined,  and  the 
originally-inclined  field  on  the  right  side  becomes  level. 
The  effect  is  net  as  dramatic  as  before. 

The  mixed  layer  quickly  responds  to  the  applied  heating. 
By  12  hours  the  mixed  layer  has  re-formed  in  the  far  field 
at  40  m.  The  shallow  section  within  the  front  is  advected 
with  the  plume  cf  the  front  and  moves  into  the  cooler  (left) 
side.  The  original  shallow  section  of  the  mixed  layer  depth 
at  the  right  boundary  of  the  front  has  deepened  to  the  mean 
far- field  value  prior  to  hear  12,  "filling  in"  at  the  rear 
as  the  shallow  frontal  peak  is  advected.  This  shallow  "bub¬ 
ble"  progresses  to  the  left  with  the  front  causing  shallow¬ 
ing  in  advance  of  the  frott  and  deepening  behind.  The  mixed 
layer  responds  with  changes  of  up  to  10  to  15  a  on  a  time 
scale  of  something  less  than  12  hours  (compare  Fig.  33  to 
Pig.  5). 

The  along-frent  velocity  maintains  a  constant  shape  and 
configuration  over  the  48-hour  integration.  The  maximum 
speed  decreases  slowly  over  time  from  100  ca/s  at  the  start 
to  6C  ca/s  by  hour  36.  The  horizontal  gradient  on  the  right 
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Figure  33. 


Front  2  Case  III  36-Hour  Hired  Layer  Depth 


side  decreases  rapidly  but  on  the  left  sida  is  prassrved 
(Fig.  34)  This  is  due  to  the  advance  of  the  plutoe  and  the 
resistance  to  nixing  that  the  heating  has  established. 

at  24  and  36  hours,  the  across-front  velocity  is  neg¬ 
ative,  with  the  maximum  cf  -10  ca/s  occurring  within  the 
aoving  pluae  cf  cverriding  wars  wafer  (Fig.  35) .  &  wave 
oscillation  in  the  v  field  is  evident,  and  has  the  same 
direction  and  time  scale  as  these  of  Cases  I  and  II.  At  12 
hours,  v  is  positive,  which  is  probably  a  response  to  the 
ageostrophic  u-coaponent  created  as  the  front  undergoes 
adjustment.  At  48  hours  the  v-field  has  a  small  positive 
and  may  be  inertial  action. 

B.  SOMMABI  —  CASE  III 

The  application  of  heating  is  able  to  generate  and  pre¬ 
serve  a  warm  pluae  which  overrides  the  cooler  water  in 
response  to  net  Ekman  transport.  Some  downwelling  is 
implied  at  the  leading  edge  of  the  front.  The  isotherms 
show  a  "kirk"  and  an  ageostrophic  u-velocity  is  there.  Seme 
of  the  surface-warmed  water  is  mixed  into  the  ambient  pool. 
On  the  right  "upstream"  side,  upwelling  occurs  as  cooler 
water  is  brought  up  tc  replace  that  lost  by  transport. 
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The  mixed  layer  depth  shallows  to  a  constant  depth  outside 
of  the  frontal  area  under  the  uniform  heating,  and  the  shal¬ 
lower  peak  existent  within  the  front  is  advected  with  the 
front  in  a  response  to  the  nass  transport.  The  across-front 
velocity  generally  shows  a  maximum  at  the  surface  within  the 
front  the  direction  of  net  Ekaan  transport,  though  reversals 
in  direction  cccur  with  tiae  as  a  result  of  ageostrophic 
responses . 
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X.  MODEL  RESULTS  FOR  PROMT  g  gV 


A.  ANALYSIS  AND  DISCUSSION 

The  combined  effects  cf  surface  heating  and  a  negative 
wind  stress  act  to  push  the  upper  150  m  of  the  front  intc 
the  warmer  water  on  the  right  side.  Vertical  mixing  is  pro¬ 
nounced  as  seen  by  the  vertical  orientation  of  the  surfaced 
isotherms  (Fig.  36)  .  This  result  is  expected  as  convective 
mixing  is  ccupled  with  wind  stirring.  Diffusion  and  hori¬ 
zontal  mixing  cause  the  frontal  horizontal  temperature  gra¬ 
dient  to  decrease  continually  over  the  48-hour  integration. 
Here  too,  the  entire  field  develops  a  dowrslope  to  the  right 
as  in  the  preceding  three  cases. 

The  mixed  layer  is  initially  shallowed  to  around  40  m  by 
the  applied  heating,  and  the  shallow  section  located  within 
the  front  breaks  down  and  approaches  the  level  of  the  over¬ 
all  mixed  layer  (Fig.  37) .  This  result  is  expected  based  on 
previous  cases  which  combined  the  action  of  convective  mix¬ 
ing  and  wind  stirring. 

The  u-velccity  field  has  an  ageostrophic  component  at 
the  left  boundary  as  a  response  to  the  adjustment. 
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Fiaure  36.  Front  2  Case  11  36-Hour  Temperature.  Contour 
""  interval  is  0-3  deg  C. 
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H  at  hour  36 


Figure  37.  Frcnt  2  case  IV  36-Hour  Mixed  Layer  Depth 


Basically  the  u-velocity  decreases  it3  maximum  speed  and 
weakens  its  bcrizcntal  gradient  over  the  48-hoar 
integration.  The  across-front  v-velocity  is  in  the  direc¬ 
tion  cf  transport  at  12  and  48  hours  and  extends  free  the 
surface,  through  the  mixed  layer  depth,  and  down  to  a  depth 
below  100  s.  Once  again,  cross-frontal  mixing  occurs.  At 
24  heurs,  the  v-field  is  reversed  over  the  entire  upper  2C0 
a.  At  36  hours,  the  v-field  exists  only  in  the  shallow  sur¬ 
face  layer,  and,  by  48  hours,  the  field  returns  to  the  posi¬ 
tive  direction.  Some  type  of  oscillation  is  evident  in  the 
v  field  (refer  to  Pigs.  3e,  39,  40,  and  41).  At  hour  12, 
the  field  is  positive  everywhere.  At  hour  24,  it  is  neg¬ 
ative  everywhere.  At  hour  36,  it  is  negative  in  the  shallow 
depths  only.  At  48  hours,  it  is  again  positive  everywhere. 
This  pattern  cccurred  in  the  previous  three  cases,  and  it  is 
stronqly  suggestive  of  seme  type  of  oscillation  in  the 
v-coaponent  both  in  direction  and  in  depth-extent. 

B.  SOBBABY  --  CASE  I? 

The  adjustment  of  the  front  under  a  negative  wind  stress 
and  surface  heating  is  affected  by  both  convective  mixing 
and  wind-stirring  in  the  vertical.  Cross-f rontal  mixing  is 
evident  in  the  horizontal.  The  front  is  "pushed  back"  into 
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the  warier  water  on  the  right  by  the  transport  of  cooler 
water.  The  front  diffuses  over  tine. 

The  wired  layer  depth  shallows  outside  of  the  front  due  to 
the  surface  heating,  and  the  shallow  portion  within  the 
frcnt  is  ercded  through  airing  and  approaches  the  depth  of 
the  far  field.  The  u-velccity  follows  an  uneventful  pattern 
of  diffusicn  and  a  decrease  in  marinua  speed  as  the  front 
adjusts.  The  direction  of  the  v-velocity  correlates  to  the 
direction  cf  net  Ekman  transport  only  at  12  and  48  hours. 
Between  those  times,  an  oscillatory  pattern  in  direction  and 
depth  is  observed.  This  is  the  same  pattern  which  was 
observed  it  the  v-fields  cf  Cases  I,  II,  and  III.  The 
v- field  indicates  cross-frontal  airing  as  it  ertends  over  a 
great  depth.  It  is  independent  cf  the  mired  layer  depth. 
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V  at  hour  24 


Pigure  39.  Frcnt  2  Case  IV  24-Hcur  V-Velocity. 
interval  is  5  ca/s. 


Contour 


Depth  (m) 


V  at  hour  48 


Piqure  41.  Front  2  Case  IV  48-Hour  V-Ve locity.  Contour 
interval  is  5  cm/s. 


XI.  COSCIOSIC NS  AND  RECCMflSND ATIONS 

A.  GENERAL 

This  thesis  has  examined  the  unsteady  response  of  two 
oceanic  fronts  to  the  local  ataospheric  forcings  of  wind 
stress  and  surface  heat  flux.  The  sinplicity  of  Front  1 
enabled  identification  of  basic  processes  stiaulatad  by  the 
forcings,  processes  which  were  not  so  clearly  seen  in  Frcnt 
2.  Front  2  was  a  acre  realistic  simulation  of  an  oceanic 
front,  and,  hopafully,  it  will  serve  as  a  stepping  stone  for 
further  investigations  into  fronts  and  testing  of  this  fron¬ 
tal  model. 

There  are  a  few  very  general  properties  that  exist  for 
both  Fronts  1  and  2,  and  which  bear  out  our  hypotheses 
stated  in  Chapter  2.B.: 

1)  Hind  direction  plays  a  primary  role  in  frontal 
adjustment.  In  the  cases  where  warmer,  lass  dense  water  is 
transported  tc  cooler,  denser  water,  the  structure  of  the 
front  is  preserved.  The  mixed  layer  depth  is  also  pre¬ 
served,  and  little  cross- frontal  mixing  occurs. 

2)  In  the  cases  where  cooler,  denser  water  is  trans¬ 
ported  to  warmer,  less  dense  water,  the  horizontal  tempera¬ 
ture  gradient  is  diffused,  the  mixed  layer  depth  within  tbs 
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front  breaks  down  and  deepens,  and  such  cross-frontal  mixing 
occurs,  ccnvective  mixing  occurs  along  with  wind  stirring. 

2)  &  uniform  constant  surface  heat  flux  shallows  the 

mixed  layer  outside  of  the  frontal  zone  to  a  depth  which  is 
propcrtional  tc  the  amount  of  heat  flux.  The  combination  of 
a  +1  dyne/cm^  wind  stress  and  a  surface  heat  flux  of  -0.004 
deg-cn/s  nixes  heat  down  to  a  depth  of  the  initial  mixed 
layer  and  leads  tc  stratification  there.  Stratification  of 
the  isotherms  at  the  surface  did  not  occur.  For  the  cases 
in  which  tie  heat  flux  is  applied,  turbulent  mixing  is 
restrained  in  the  vertical,  which  agrees  with  the  results  of 
Niiler  (1975)  and  Mellor  and  Durbin  (1975). 

4)  In  areas  of  the  frontal  zone  where  ageostrophic 
velocities  are  created,  vertical  mixing  is  increased. 

5)  Mixed  layers  deepen  rapidly  within  zhs  first  12 
hours  under  wind  stress  alone  as  the  fronts  undergo  initial 
adjustment.  Beyond  12  hours,  adjustments  are  confined  to  an 
area  around  tie  the  basic  depth  attained  during  the  initial 
adjustment  period.  Miller  (1975),  Mellor  and  Durbin  (1975) 
and  Oe  Szceke  (1980)  all  reported  this  phenomenon. 

6)  A  time-dependent  pattern  in  the  direction  and 
depth-influence  of  the  across-front  velocity  exists.  it  is 
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probably  due  tc  inertial  oscillations  created  by  the  atmos¬ 
pheric  forcing. 

It  was  surprising  to  see  the  downward  bulge  in  the  mixed 
layer  at  the  cold-side  boundary  of  both  fronts  in  the  wind 
stress  only  cases.  It  is  unclear  as  to  why  this  happens, 
though  soae  speculative  reasons  vers  provided. 

In  Front  1,  the  formation  of  cores  of  positive  v-veicc- 
ity  which  are  symmetric  about  the  mixed  layer  depth  are 
intriguing.  Inertial  oscillations  may  provide  an 
explanation. 

The  tilting  down  to  the  right  of  the  entire  temperature 
field  in  Front  2  in  all  four  cases  is  a  major  adjustment, 
and  one  which  is  independent  of  wind  stress  direction  and 
heating.  It  would  appear  that  this  tilt  is  a  response  to 
the  re-orientation  of  the  oceanic  current  fields. 

B.  BECOB  HENC ATICKS 

There  is  both  further  investigation  that  can  be  done 
with  the  results  of  this  thesis,  as  well  as  new  work  which 
is  suggested  by  these  results.  The  results  of  these  eight 
model  runs  could  bm  further  quantified  by  using  a  finer  con¬ 
tour  interval  for  the  fields,  and  by  producing  fields  for 
gradients,  agecstrophic  velocities,  and  the  like.  The 
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inertial  actions  coaid  be  investigated  to  encompass  magni¬ 
tudes,  directions,  depths  and  areas  affected.  It  would  be 
instructive  to  examine  <y,z)  streamf unctions.  The  model 
could  be  tested  by  conducting  parametric  sensitivity  studies 
of  the  constants  used  in  the  governing  equations  and  mixed 
layer  equations. 

New  work  might  include  using  different  combinations  of 
these  atmospheric  forcings  both  in  magnitude  and  direction. 
Inclusion  of  a  wind  stress  in  the  y-dir9Ction  would  cer¬ 
tainly  prove  to  he  of  interest,  as  would  a  diurnal  heating/ 
cooling  cycle,  ether  atmospheric  forcings  could  be  applied. 
Evaporation  and  precipitation  might  be  incorporated.  The 
equation  of  state  could  be  expanded  to  Include  salinity,  and 
an  initial  salinity  field  could  be  inserted  and  monitored 
over  time.  Cf  course,  a  three-dimensional  version  of  this 
model  would  be  valuable.  More  simulations  of  observed 
fronts  using  data  acquired  in  the  field  would  provide  fur¬ 
ther  insight  into  the  applicability  of  the  model.  Data 
acquired  by  Jchannessen,  et  al  (1977)  of  the  Haltese  front 
would  be  a  recommended  starting  point. 

Computer  central  processor  unit  (CPU)  time  was  extensive 
for  each  of  these  model  runs.  They  each  averaged  over  four 
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;  hours  of  CFQ  tins  to  run,  providing  no  stopping  and 

i 

ing  were  required.  Such  a  nodal  is  far  from  bring  useful 
operationally  fcr  real-time  assessment  of  ocean  structure  to 
provide,  for  example,  input  for  acoustic  forecasting  acdels. 
However,  fcr  basic  research,  this  model  appears  to  be  a  pow¬ 
erful  tool  for  deciphering  the  complexities  of  ocean  frontal 

iynanics  and  thermodynamics. 
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